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ABSTRACT 

An analysis has been made of some aspects of atmos- 
pheric sensing by laser radar spectroscopy. A discus- 
sion of the laser radar equation is presented which leads 
to methods for obtaining atmospheric information by 
spectroscopic means. Calculations are presented which 
investigate the potential of Doppler, Raman and differ- 
ential absorption techniques. It is concluded that at- 
mospheric measurements based upon Doppler and Raman 
approaches are restricted to qualitative studies at the 
present time primarily because of the limitations im- 
posed by existing laser sources and detectors. It is 
probable that these limitations will be removed in the 
near future and it is recommended that exploratory 
experiments be continued. The differential absorption 
method for obtaining the spatial distribution of atmos- 
pheric gases is well developed for such gases as water 
vapor and oxygen where the basic laser frequency lies 
close to the center of the particular absorption line. 
Further investigation is recommended in the field of 
laser frequency shifting. Additional information is 
also required on the shape of the absorption line, 
particularly in the region where the line reflects both 
pressure and Doppler broadening. 


1. INTRODUCTION 

The purpose of this paper is to consider the potential of the laser radar as a 
ground based remote sensor of the earth's atmosphere. Perhaps the simplest ap- 
proach to this subject may be had by a consideration of the transfer equation which 
relates the power transmitted by the laser to that which is captured by coaxial re- 
ceiving optics. 
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Radiance Scatter cross Scattering Transmis- Capture 
at Z section per volume sion of at- solid 
unit length mosphere angle 

= received power 
P T = transmitted power 
A^ = receiver aperture 
Z = slant range 

(Tj. = molecular backscatter cross section 
n(Z) = molecular number density 
AZ .= pulse length 


T = two way atmospheric transmission 

0 = receiver beam width 


( 1 ) 


When this equation is applied to the atmosphere it takes the normalized form: 

z2 P^ = AZ n < Z ) + ^2 J i(r,X)N(r) drj T (2) 

o 

Z Z Z 

T = ex P ~2[ J dZ + | /3 A dz + | p(Z) dz] (3) 

o O o 

In this equation the following definitions hold: 

= average irradiance at the receiver at wavelength X 
= pulse power of the transmitter at wavelength X 
= slant range to the pulse volume 
= laser wavelength 
= Mie scattering function 
= molecular number density at range Z 
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N(r, Z) 
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spectral distribution of aerosol particles at range Z 
molecular scatter cross section/unit volume 
aerosol scatter cross section/unit volume 
molecular absorption cross section (area/unit mass) 
277 molecular polarizability 
atmospheric density 


la. Atmospheric Constituents 

The atmosphere is composed primarily of nitrogen (78. 5%) and oxygen (20. 95%). 
These gases are well mixed and maintain the same percentages on a volume basis 
throughout the troposphere and stratosphere. The major variable gases of meteor- 
ological interest are water vapor, ozone, and carbon dioxide. These gases can vary 
rapidly both in time and in space. 

The surface air contains large numbers of small particles whose sizes range from 
10“ 2 microns to 2 microns. Total number densities may be as large as 10^ per cubic 
centimeter. The bulk of these particles have radii less than 0. 1 micron and a typical 
count for particles greater than 0. 5 micron would be 1 per cubic centimeter. The par- 
ticle density normally decreases with height in an exponential manner as shown in 
Figure 1. The refractive index of these particles is complex with a real part in the 
range 1.5 to 1.7. little is known of the spectral variability of either the real or 
imaginary portions of this quantity. 



Figure 1. Variation of particle concentration with height. 
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lb. Atmospheric Spectral Properties 

Elterman has computed the sea level particle scattering coefficient as a function 
of wavelength. These values have been plotted by Breece, et al. (1966) and are pre- 
sented here as Figure 2. It should be noted that the scattering coefficient changes 
very slowly in the region below 1 micron. At wavelengths near 10 microns the aerosol 
particles will be in the Rayleigh range and the wavelength dependence will be A. -4 . 



Figure 2. Variation of aerosol particles scattering coefficient with wavelength. 


The magnitude of the particle backscatter at A = . 7p is an order of magnitude 
greater than that due to molecular scatter at ground level. However, due to the rapid 
decrease of particle density with height the two terms become essentially equal at 3 
kilometers and remain so throughout the troposphere as shown in Figure 3 (Breece, 
et al . 1966). 

The strong absorption bands of water vapor, carbon dioxide and ozone in the 
infrared and of ozone in the ultraviolet regions of the spectrum are well known to the 
meteorologist as they play important roles in the thermodynamics of the atmosphere. 

It is normally assumed that the atmosphere is essentially transparent from 0. 3 to 0. 9 
microns. However, reference to high resolution telluric spectra shows that narrow, 
intense absorption lines due to water vapor, oxygen, ozone are dispersed throughout 
this region of the spectrum. A sample of such a spectrum taken at N. Y.U. by Bradley 
is shown in Figure 4. It should be noted that the atmospheric absorption for water 
vapor (one way)at 6943. 80A is nearly 60% and the line width is essentially 0. 05X at half 
width. Reference will be made to this isolated line in a latter portion of this paper. 
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Figure 3. Ratio of particle to Rayleigh scattering coefficients as a function of height. 


lc. Energy Background Level , 

In the discussion which follows it will be necessary to compute signal to noise 
ratios of various proposed systems in order to establish their feasibility. Generally, 
these systems, if they are to be considered useful from a meteorological standpoint, 
must be operable during day as well as night hours. Consequently, the noise contri- 
bution due to the presence of background radiation must be known. The background 
radiation may be, computed from the relationship 

P b (M =A r J( X ) A CO 

where 

Pg(A) = background radiation (watts) 

O 

= receiving aperture (cm ) 

J(A) = radiance (watts /cm 2 /sr/micron) 

Ate = solid viewing angle 
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Figure 4. Telluric absorption spectrum (6933-6945X) April 10, 1967 
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A graph of J(A) versus A of blue sky observed by Vanderhei and Taylor (1960) is 
given in Figure 5. A value of 10 - 3w/cm 2 /sr/micron is characteristic of the ruby laser 
region. 


2. CAPABILITIES OF LASER SYSTEMS 

One of the goals of a remote sensing device would be to provide information equi- 
valent to that obtained by the standard radiosonde under all weather conditions. It is 
clear that the laser radar will not penetrate thick cloud decks. This limitation is im- 
posed by the available energy in laser pulses. For example, assume a cloud deck of 
300 meters thickness and of 3 gram/meter^ liquid water content. In order to penetrate 
this cloud and provide a path for the return scattered energy it would be necessary to 
evaporate the water in the path of the beam. If the beam were 0. 1 meter x 0. 1 meter 
in cross section this would require a pulse containing approximately 21, 000 joules of 
energy. Present day systems provide energy levels (Q switch) of a maximum of 10 
joules. It does not seem possible in the foreseeable future to provide energy of the 

requisite amount through combinations of increased pulse energy and repetition rates. 

\ 

What can the laser radar do aside from deriving qualitative indications of n(Z) and 
N(Z) ? We would like to operate the laser radar in a manner which would permit the 
separation of the quantities n(Z), p(Z) and T(Z). The following approaches to such 
measurements will be considered in the remainder of the paper 

n(Z) and N(Z) via Raman scattering 

p(Z) via differential absorption of scattered energy 

T(Z) via Doppler power spectrum of the scattered energy. 

v(Z) via Doppler power spectrum of the scattered energy. 


2a. Doppler Power Spectrum 

The transfer equation stated in the introduction of this paper relates the energy 
returned to the laser radar to that which is transmitted. The scattering process is an 
interaction mechanism between the traveling laser energy and the molecules and par- 
ticles which comprise the atmosphere. These constituents are in random motion and 
their velocities are governed by the Maxwell-Boltzmann relationship 


dN(v) 


1_ dv 
N \fW 


CT 

V 


exp 


where 


(v - v) 2 

2a 2 


(4) 


v = mean wind 

v = random particle velocity 

— standard deviation of the thermal wind due to Brownian motion 
N = number of particles/unit volume 

dN(v) = number of particles in the range dv at the velocity v 


186 



R. M. SCHOTLAND 


The scattered radiation although initially nearly monochromatic will be Doppler 
broadened by virtue of the velocities of the scattering components. The Doppler shift 
may be expressed for these particles as 



a a 
(r + t) 


where 

= particle velocity 
= unit vector to the receiver 
= unit vector to the transmitter 
for a coaxial system 



v. 

1 

A 

r 


$ 


This relationship may be used with the Maxwell-Boltzmann expression to trans- 
form the normalized transfer equation into a power spectrum for the backscattered 
radiation. Use is made at this point of the relationship 


cr 


2 


v 


kT 

m 


( 5 ) 


where 


k = Boltzmann's constant 
T = temperature 
m = scattering mass 


( 6 ) 


It is assumed 


simply as 
This expres- 
sion does not deviate appreciably at halfwidth from the Boltzmann form for molecular 
scatter, but does so markedly for particles larger than 0. Ip (Schotland, et al . 1967). 
The Boltzmann form also has been used by Zirkel (1966) and Breece, et al . (1966). It 
is important that this point be clarified. 


_ d H R 

P T ^ f d 


ZJ?L AZ n( Z ) 

A 2 [277R*T] 1//2 


A 2 ft 
o d 

* 

8R T 


This relationship (Schotland, et al . 1962) is plotted in Figures 6 and 7. 
that the plot has been centered about _ v 

*d " A ’ 


It should be said at this point that the spectrum may not transform 
shown. Dicke (1953) has suggested a collision narrowed Lorentz form. 
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Figure 6. Representation of molecular and aerosol power spectra. 


The normalized power spectrum for particles in the Rayleigh range is 


R 


2 ^ H R 

df. 


— o K Az V 2 N(Z) 


V2770- 


exp 


A 2 f 2 
A f d 

8a 2 


where V - volume of scattering particle. 


„ _ m - 1 

K “ m + 2 ’ m = refractive index. 


(7) 


This form of the equation is useful at the C0 2 laser wavelength (10 . 6/i). 

In order to observe the power spectrum it is necessary to limit the receiving 
aperture to the coherence area of the sources. This area, A c is of the order of A 2 /Cl 
where 1 1 is the solid angle of the scattering volume viewed from the receiver. This 
restriction implies that the beamwidths of the transmitted signal must be extremely 
narrow if practical values of A c are to be obtained. Fried (see Breece, etal. 1966) 
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Figure 7. Normalized molecular Doppler spectrum. 
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has demonstrated that the transmitting optics should be diffraction limited and equal 
in area to that of the receiver for optimum results. 


Another problem area in the application of Doppler techniques to the atmosphere 
results from the perturbation of the wavefront of the radiation by the atmosphere. The 
diameter of the receiving telescope should not exceed a value r 0 given by Fried (1965) 
as 


r 


o 


1.2 x 


10“ 3 * * * * 8 A 6/5 R- 3/5 


-6/5 

N 


(8) 


where 

A = wavelength 

R = range (two way) 

2 

C.. T = Tatarski's structure constant 


The strong dependency on A should be noted. As an example r Q for day operation 
is of the order of 7 cm for A = . 7p while it may be as large as 10^ cm at lOp. 

In order to discuss the potential of the laser radar for Doppler application, the line 
purity of typical lasers must be considered. At the present time two classes of lasers 
can be considered useful. They are the pulsed crystal and the C0 2 gas laser. 

According to Hercher (private communication) the line width of pulsed crystal 
lasers can be of the order of 30 x 10® Hz, when operated in the Q-switch mode. This 
line width is of the order of the pulse spectrum so that no improvement appears pos- 
sible except through the use of longer pulses. The energy available is approximately 

0. 1 joule in projected laser systems. 

From a discussion with equipment manufacturers it appears that the line width of a 
CC>2 gas laser is of the order of 0. 008X (2.4 x 10® Hz) at a power level of 100 watts. 

Modulation schemes have been proposed whereby the Doppler shift occurs at the 
modulation frequency. This technique has the twin advantages of strong spectral purity 
and a coherence area dictated by the modulation frequency. However, present modu- 
lation schemes are limited to 10 9 Hz (A = 30 cm). Therefore the Doppler shifts are 
generally small: f d /v = 6 Hz/m/sec. 


3. APPLICATIONS OF THE DOPPLER SYSTEM 

It has been suggested that the Doppler technique be used to measure the following 

parameters. 

1. Temperature from the width of the molecular spectrum 

2. Wind velocity via the mean Doppler shift using heterodyne techniques 

3. Velocity gradient using square -law detection of the backscattered radiation 

from two laser pulses separated in range 
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Problems 1 and 2 were analyzed by Schotland, et al . (1962) and although times 
have improved and a better understanding of the problem is available, the conclusion 
is still pessimistic. 

Consider the temperature problem. 

I = 

1(f ) = 3.4 x 10 *” for A = 0.7p 



Then at the half power point which is sufficiently removed from the aerosol contribution, 

'(I)- 1 - 


,7 x 10 


AP„ = 1. 7 x 10~ 16 x —77- A x Af 
R ,-,2 c 


a range 


Phase limitations limit A c to approximately 50 cm 2 . Assume Af = 10? 
nge of 1 km. 


Hz and 


km. 

J 

R 


1 /? n IQ 7 

AP p =1.7x10 x 10 x 10 x 50 x 10 

(Af) (Z 2 ) (A) (P) 


<y 


AP d = 8.5 x 10 11 watts 

ri 


Substitution of this power into the S/N relation assuming Pg = 0 yields S/N = 3 
which is not very satisfactory. It should be noted that with coherent detection the S/ 
-ratio is hasinallv ormivalont to that nhtainpH with Pt-* = 0. 


li CXDOUXXXXXXg jrg — \j y icius o/ in — o 

wmcn is not very satisiactory. it snouia oe notea that with coherent detection the S/N 
ratio is basically equivalent to that obtained with Pg = 0. 

An estimate of the desired signal/noise value for thermal measurement can be 
obtained by logarithmically differentiating the expression for I. 


Thus for 


X 2 f 2 
o d 

8R*T 


dl 

I 


dT 

T 



L8R T 


1 

2 


1, a 2° temperature precision at the level of 300°K requires a pre- 


cision in the measurement of I of 1/300. This level of accuracy could be obtained at 
a maximum range of only 10 m. 


Bistatic systems using c. w. lasers such as the Argon or the CO2 units can be used 
with synchronous detection techniques for the temperature measurements. The dif- 
ficulty with this approach is that with present capabilities the integration times neces- 
sary to yield reasonable S/N ratio are of the order of minutes. Another method which 
can be used to obtain the change in Doppler width involves the direct use of an inter- 
ferometer. The advantage of this approach is that the conditions on phase coherence 


\ 
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are relaxed - permitting much larger receiving apertures. The interferometer pro- 
vides information only on the spectral width and therefore this method cannot be used 
for velocity studies. 

Doppler techniques to measure the mean wind and the radial wind gradient are 
basically dependent upon the ability of the system to "see" these shifts in the presence 
of the spectrum variability of the laser source. The Doppler shift sensitivity may be 
expressed as f^/v = 2/X. 

X . 7/i 10. em 


Y v 


3 x 10^m ^ 


2 x 10 5 m -1 


If a velocity measurement accurate to 1 m/sec is desired, the expected shift will 
be of the order of 0. 1 of the spectral noise of state of the art laser systems. This un- 
certainty does not lead to a useful measurement. At the present time a laser is not 
available which will deliver 0. 1 joule with a bandwidth of 30 x 10 6 Hz. It appears that 
such a system could be developed with existing technology. 

Recently Fiocco and DeWold (private communication) performed an experiment in 
which they used a pressure scanned Fabry-Perot interferometer to observe the 
Doppler spectrum of an aerosol. Although this system cannot be used as it stands for 
atmospheric probing, the work does demonstrate the reality of the Doppler spectra. 

Similar conclusions were also reached by Cummins, et al . (1964) who used optical 
mixing techniques to study the Doppler power spectrum in liquids. 

It is recommended that an experimental program be initiated to study the thermal 
Doppler spectrum. In particular the shape of the spectrum should be considered in 
relationship to the Lorentz versus Gaussian forms. 


3b. Raman Scattering Technique 

In general, the radiation which is backscattered is returned at a wavelength 
similar, except for Doppler shift, to the existing wavelength. At times the scattering 
molecule will undergo a change of state during the scattering process and the scattered 
energy will be returned at a different wavelength. This process is known as Raman 
scattering. The average Raman cross section is of the order of 4 x 10 -3 that of the 
Rayleigh cross section for atmospheric gases. The changes of molecular state are 
related to the vibration and rotation spectra of the molecules and the magnitude and 
wavelength shift will be specific for a given molecule. Since the Raman cross section 
(non-stimulated) is small, the technique will probably be useful only for the major 
atmospheric components such as N 2 and 0 2 . 

An example of the signal to noise ratio in the measurement of the vertical distri- 
bution of C>2 is calculated below. It is assumed that a ruby laser operating at 6943& 
has been frequency doubled with an efficiency of 25 percent. The V = 0 to v = 1 vibra- 
tion state corresponding to a line shift of 1550 cm - -*- will be used in the computation. 
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System parameters 

A_ = 4 x 10 cm 

XV 

All = 10 ^ sr 
AA = 4 x 10 ^ ,u 

= 2 x 10 ^ watts /sr/jLt/cm^ 

P T = 0. 5 x 10 8 watts 
Af = 10 6 hz 

C = photomultiplier responsivity = . 15 amp/w. 

The signal to noise ratio for shot noise limited system such as a photomultiplier is 
given by the following expression: 

S _ P R C T 


N r „ -1/9 W 

[2e Af CT(P r + P B )f /Z 

where e = electronic charge 

Substitution of the above parameters into the S/N relationship results in the following 
values: for ranges Z = 3 and 6 Km. , S/N a 80 and 18, respectively. 

It appears that the measurement is feasible. Some work has been done by Cooney 
at R.C.A. (private communication) using a straight ruby laser system and filter suit- 
able for the N 2 molecule. His results have not been published but it appears that he 
has observed the Raman profile. 

There is a difficulty with this model. While Raman scattering is specific for a 
particular molecule, the attenuation function is not. Thus the returned signal will be 
modulated by the total atmospheric cross section. A possible solution to this problem 
would be to monitor the unshifted backs catter (or an additional line) and use this quan- 
tity as a correction to the Raman measurement. 

Consider the relationship 


L r + 

P T 2 [RAY 


for the unshifted line 


RAM 


- 7T ^ ct ram J 7-2 T 2 - t for the Raman line 


T 1 CT P 

= k 1 + k 0 — — 

l 2 1 2 CT RAY 


i. e. the ratio of particle to 
Rayleigh cross section is 
available 
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Cooney (private communication) points out that temperature information is incor- 
porated in the ratio of the intensities of the Stokes and anti -Stokes lines. It is not 
clear at the present time if the intensity of the anti -Stokes line is of sufficient magni- 
tude to insure a useful measurement. 


3c. Differential Absorption of Scattered Energy 


The expression for atmospheric transmission which appears in the radar transfer 


equation includes a term exp 



pK^ dz . 


This function represents the molecular 


o 

absorption due to atmospheric gases. The atmospheric spectrum presents an abun- 
dance of narrow lines (O.lA) of moderate strength due to H 2 0, 0 2 , O3, and CC^. 
Figure 4 represents a portion of the telluric spectrum from 6934 to 6944X obtained by 
Bradley. The predominant lines in this region are due to H 2 0 and 0 2 . A technique 
has been developed at N. Y. U. over the past 10 years to measure the vertical distri- 
bution of absorbing gases by measuring the differential absorption of the backscattered 
radiation from a pulse emitted by an optical radar. The technique was initially applied 
to O3 in 1958 using a pulsed Xenon source. However, the results were marginal due 
to the poor spectral brightness of the source. The development of the laser has pro- 
vided an energy source which is ideal for this type of measurement. Basically the 
laser is tuned to an absorption line center and then to an adjacent window region. 


An example of simulated radar data is presented in Figure 8 which represents 
the transfer equation evaluated for 6943. 80A water vapor line. Line parameters are 
due to Benedict. 


The inversion of the optical return to obtain the vertical profile of water vapor 
proceeds as follows: 


R(Z) 


P R(* l ) F T(X l ) 

Pr <V Ft <V 


exp 


- 2 / p y ( Z)K X (X.) dz 


5z ln R(Z) 


-2 P v (Z)K w (X.) dz 


••• P V < Z > = 


2kV> W R(Z) 


w' 1' 


(13) 

(14) 

(15) 


The accuracy of the determination of the water vapor profile by the differential 
absorption technique is limited primarily to the uncertainties in the measurement of 
the gradient of the power return and by the incomplete knowledge of the water vapor 
absorption coefficient. The power measurement, at least from the first four kilo- 
meters of the atmosphere, does not represent a problem area. 


The variance of the deduced water vapor density may be expressed in terms of 
the variances of the absorption coefficient and returned power gradients as follows: 


2 

w 


- P 2 (Z) 
t ,2 



a 


2 

(dP L /dZ) 



( dp w /az) 


(16) 
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where 



= water vapor variance 


- the absorption coefficient variance 


cr 


2 

0P/dZ) = 


the power gradient variance 


The variance of the absorption coefficient is due to two factors. The first factor 
is that an uncertainty exists in the measurements of the line parameters and in the 
effects of pressure and temperature broadening upon line shape. The values for the 
line width and strength according to Benedict are thought to be correct to ±33% while 
the line strength has an uncertainty of ±25%. The second factor that contributes to the 
variance in the absorption coefficient is the uncertainty in the frequency of the laser. 
From the work that has been done at N. Y. U. , it appears that a laser frequency setting 
uncertainty of ±0. 01A is possible. If the laser is tuned to line center and the half 
width of the water vapor is taken to be 0. 04&, then the fractional uncertainty in Kx is 
of the order of 5 percent. However , if the laser should be tuned to the inflection point 
of the curve, the fractional uncertainty can be as large as 29%. 

The variance of the water vapor density due to the power gradient term arises 
primarily from the shot noise produced by the photomultiplier-detecting system. The 
variance of radiant noise power may be written 

< = f ( P r + P b) A£ < 17 > 

where 

e 

C 

Af 


The variance 


2 ~ N 

CT (dp/dP) Az 2 ( 18 ) 


= electron charge 
= photomultiplier responsivity 
= measurement system bandwidth 
= background power 

dP 

of can be expressed as 


Substitution of the expression for variance of (&P/dZ) into the equation for the 
variance of p w yields 


ct 2 *2 + g_ Af 

w A K _ * ..2 . 


AK 


rak 2 AZ 2 


P L + P B 
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P + P_ 
w B 


P 2 
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(19) 
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The value of 0^, has been computed using the model atmosphere of Figure 8 and the 
following constants : 

Pg = 5 x 10 ® watts 

7 -9 

Pj. = 2 x 10 watts (over 25 x 10 sec) 

A = 10 3 cm 2 

Af = 2 Az/c = 4 

Az = 120 meters 


The results are plotted in Figure 9 separately for the contributions due to an d 

to or . The uncertainty in cr was assumed to arise solely from the frequency stability 
of the laser. It is apparent that in the first two kilometers the uncertainty in O' arises 
from cr^.. However, the contribution from increases rapidly above two kilo- 

meters and is the dominating factor at three kilometers. The uncertainty in p w due to 
°dP/dZ can re( iuced by increasing both the output energy of the laser and the aperture 
of the receiver. 


The major problem area in developing a field laser radar system of this type cen- 
ters about the frequency stability of the laser. It is necessary to center the laser on 
line center with an accuracy of ± . OlX. Initial experiments were undertaken using the 
thermal tuning characteristic of the ruby laser to locate the line center. The thermal 
sensitivity of a ruby laser is 0. 065X/°C which implies a rod thermal stability of 0. 1°C 
if the wavelength tolerance of . OlX is to be maintained. This proved difficult to achieve 
at high laser repetition rates and consequently resonant cavity frequency stabilizers 
are now used. These devices also permit rapid frequency shifting of the laser. A 
sample of R(Z) data for water vapor obtained using a thermally tuned laser is given in 
Figure 10. 

It is possible in principle to generate both the window and line center wavelength 
by using an appropriate Raman shifter as suggested by Dobbins and La Grone, (1967). 
However, no calculations are available in which the power output and line width of the 
appropriate Raman shifted energy are available. This work should be continued since 
Raman shifting will provide access to a variety of atmospheric lines. 
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